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Abstract 

We consider a supergravity (SUGRA) scenario, with universal scalar and gaugino 
masses at high scale, with a right-chiral neutrino superfield included in the spectrum. 
Such a scenario can have a lightest supersymmetric particle (LSP) dominated by the 
right sneutrino and a stau as the next-to lightest supersymmetric particle (NLSP). Since 
decays of all particles into the LSP are suppressed by the neutrino Yukawa coupling, 
the signal of supersymmetry consists in charged tracks of stable particles in the muon 
chamber. We demonstrate how a neutralino decaying into a tau and the stau-NLSP can 
be fully reconstructed over substantial areas in the SUGRA parameter space. We also 
suggest event selection criteria for eliminating backgrounds, including combinatorial 
ones, and use a new method for the extraction of the mass of the stau-NLSP, using its 
three-momentum as obtained from the curvature of the charged track. 
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1 Introduction 



Searches for supersymmetry (SUSY) [1, 2, 3] at the Large Hadron Colhder (LHC) are largely 
based on signals with missing transverse energy (^t)- This is because SUSY, in its 7?.-parity 
conserving form (with 7?.-parity defined by TZ — (_)3S+l+2S'j^ offers the lightest supersym- 
metric particle (LSP) which is stable, and, if electrically uncharged and weakly interacting, 
is potentially a cold dark matter candidate. The hghtest neutralino (xi) turns out to be 
the LSP in most theoretical models. All SUSY cascades at colhder experiments should 
culminate in the pair-production of the LSP within the detector itself. The neutral and 
non-strongly interacting character of the LSP results in its invisibility at colliders, and thus 
a large energy- momentum imbalance, together with energetic jets and/or leptons emerges as 
the characteristic signal of SUSY containing a dark matter candidate. 

It should be remembered, though, that the above possibility is not unique. Apart from the 
lightest neutralino, the left-chiral sneutrinos in the minimal SUSY standard model (MSSM) 
can in principle be a dark matter candidate as well. This is, however, strongly disfavoured 
by direct dark matter search experiments, because the SU(2) interaction of a left-chiral 
siicutrino (as opposed to the U(l) or Yukawa dominated interaction of a neutralino) gives 
rise to unacceptably large cross-sections of elastic scattering with dark matter detectors. In 
addition, a left-chiral sneutrino LSP is difficult to accommodate in a scenario where the 
SUSY breaking masses evolve from 'universal' scalar and gaugino mass parameters at a high 
scale [6]. 

The situation changes if one has right-chiral neutrino supcrficlds in addition, a possibility 
that often haunts us as evidence piles up in favour of neutrino masses and mixing [4, 5]. It 
has been shown in some recent works [8] that such a right chiral sneutrino may pass off 
as a dark matter candidate without any contradiction from available data [9]. Since the 
right-chiral sneutrino has no gauge interaction, the only way it can interact with matter is 
via neutrino Yukawa coupling, the strength of its interaction is too feeble to be seen in dark 
matter search experiments. In such a case, the next-to-hghtest SUSY particle (NLSP), too, 
has an excruciatingly slow rate of decay into the LSP dominated by right-chiral sneutrino 
states. Consequently, the NLSP is stable on the scale of colhder detectors, and, in cases 
where it is a charged particle, the essence of the SUSY signal hes not in but in charged 
tracks due to massive particles, seen in the muon chambers. 

Such stable charged particles can in principle be distinguished from muons through a 
number of techniques. These include the measurement of time delay between the inner 
tracking chamber and the muon chamber, the degree of ionisation, and also more exotic 
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proposals such as the absorption of the stable particles in a chamber which can be subse- 
quently emptied underground to observe the decays [10]. While these are all of sufficient 
importance and interest, some of us have shown in earlier works [11, 12] that there are some 
very good kinematic discriminators for such stable charged particles, which make the signals 
practically background-free for both stau and stop NLSP. Event selection criteria based on 
the transverse momentum pt of the tracks, in conjuction with quantities such as the scalar 
sum of all visible transverse momenta and the invariant mass of track pairs, are found to be 
useful in this respect. In this work, we perform a detailed simulation of signals, backgrounds 
and mistags to show that the masses of neutralinos can be reconstructed to a high level of 
precision for a scenario with f NLSP and an LSP dominated by the right-chiral sneutrino 
of the third family. We use the technique of tau reconstruction for this purpose. Also, we 
depend on neither ionisation nor time delay for extracting the mass of the stable stau, but 
rather obtain it using an algorithm that depends on event-by-event information on two taus 
and two stable tracks in the final state. 

It should be mentioned that the signal discussed here as well as the reconstruction tech- 
nique advocated by us is not limited to scenarios with right-sneutrino LSP alone. One can 
have stable staus in SUSY, when, for example, one has a gravitino LSP in a supergravity 
(SUGRA) model [13]. They can be envisioned in gauge- mediated SUSY breaking theories 
as well [14]. In the MSSM, too, one can have the so-called co- annihilation region of dark 
matter, where a stau and the neutralino LSP are closely degenerate, leading to a quasi-stable 
character of the former [15]. It should be emphasized that our suggested procedure is ap- 
plicable to all of these cases. What we find as a bonus is that scenarios with stau NLSP 
and right-sneutrino LSP occur rather naturally in a SUGRA model [8, 11] with a universal 
scalar mass which is the origin of the right-sneutrino mass as well. Thus the mere addition 
of a right-sneutrino superfield, which is perhaps the most minimal input to explain neutrino 
masses, can turn a mSUGRA theory into one with a stau NLSP and a sneutrino LSP. Thus 
one can identify regions in the mo — M1/2 plane of the theory, where the reconstruction of 
unstable neutralinos is feasible at the LHC. 

In section 2, we discuss the scenario under investigation as well as the super particle 
spectrum, and motivate the choice of benchmark points used for demonstrating our claims, 
in the context of a supergravity scenario. The signal looked for, the corresponding standard 
model backgrounds and the event selection criteria chosen by us are discussed in section 3. 
Section 4 contains discussions on the various steps in reconstructing neutralinos. The regions 
in the mo — M1/2 plane in a SUGRA scenario, where neutralino reconstruction is possible in 
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our method, are also pointed out in this section. We summarise and conclude in section 5. 



2 Right sneutrino LSP in supergravity 

The supcrpotcntial of the minimal SUSY Standard Model (MSSM) [16] is given (suppressing 
family indices) by 

Wmssm = yiLHdE'^ + VdQHdD- + VuQHuU' + f^HaHu (1) 

where Ha and Hu, respectively, are the Higgs superfields that give mass respectively to the 
Ts — —1/2 and T3 = +1/2 fermions. y's are the strengths of Yukawa interactions. L and 
Q are the left-handed lepton and quark superfields respectively, whereas E'^, and U'^, in 
that order, are the right handed gauge singlet charged lepton, down-type and up-type quark 
superfields. /i is the Higgsino mass parameter. 

As has been already mentioned, the MSSM must be additionally equipped to explain 
non-vanishing neutrino masses. Phenomenologically, the simplest (though not theoretically 
the most satisfying) way to do so is to assume neutrinos to be of Dirac type and simply 
add one right-handed neutrino superfield to each family. The superpotential of the minimal 
SUSY standard model is thus extended by just one term per family, of the form 

= y.HuLi>% (2) 

However, having such small Dirac masses for the neutrinos would imply that the neutrino 
Yukawa couphngs (y,^) are quite small (~ 10~^^). The above term in the superpotential ob- 
viously implies the inclusion of right-handed sneutrinos in the (super) particle spectrum, and 
these sneutrinos will have all their interactions proportional to the corresponding neutrino 
masses. Thus the dominantly right-handed eigenstate of the tau-sneutrino might become a 
viable dark matter candidate, without coming into conflict with dark matter search limits, 
thanks to its extremely feeble strength of interaction with all matter. 

Interestingly, scenarios where the MSSM is embedded in a bigger, high-scale framework 
for SUSY breaking can support the above situation. The most commonly invoked scheme is 
minimal supergravity (mSUGRA) where all scalar (gaugino) masses at low-energy arise from 
a universal mass parameter mo{Mi/2) at a high scale where supergravity, or local SUSY, is 
broken. If one adds a right-chiral neutrino superfield, then the right-sneutrino mass may 
be assumed to originate in the same parameter mo- As we shall see below, this causes the 
physical state dominated by the right-chiral sneutrino to become the LSP. It has been shown 
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that such a possibihty is consistent with all experimental bounds [7] and also compatible 
with the dark matter density in the Universe [8, 9]. 
The neutrinos masses can be schematically shown as 



where tan/3 is the ratio of the vacuum expectation values of the two Higgs doublets that 
give masses to the up-and down-type quarks respectively. The actual mass eigenvalues will 
of course depend on the Yukawa coupling matrix. This, however, does not affect the collider 
signals of the SUSY scenario under consideration here, as the interaction strengths of the 
dominantly right-chiral states arc always very small in magnitude. 

Upon inclusion of right-chiral neutrino supcrficld into the SUGRA framework, the su- 
perparticle spectrum mimics the mSUGRA spectrum in all details except for the identity of 
the LSP. As has been already mentioned, SUSY breaking in the hidden sector at high scale 
is manifested in universal soft masses for scalars (mo) and gauginos (M1/2), together with 
the trilinear (A) and bilinear (B) SUSY breaking parameters in the scalar sector (of which 
the latter is determined by electroweak symmetry breaking conditions). Masses for squarks, 
sleptons and gauginos, all the mass parameters in the Higgs sector as well as the Higgsino 
mass parameter /i (up to a sign) are determined, once the high scale of SUSY breaking in 
the hidden sector {O ~ 10^^ GeV) are specified. Neglecting inter-family mixing, the mass 
terms for sneutrinos arising in this manner are given by 



where A^, is the term driving left-right mixing in the scalar mass matrix, and is obtained by 
running of the trilinear soft SUSY breaking term A [17]. The Yukawa couphngs can cause 
large splitting in the third-family squark and sleptons masses while the first two families are 
more closely degenerate. On the other hand, the degree of left-right mixing of sneutrinos, 
driven largely by the Yukawa couplings, is extremely small. 
The sneutrino mass-squared matrix is thus of the form 



where Mi is the soft scalar mass for the left-handed sleptons whereas the Mj>^ is that for the 
right-handed sneutrino. In general, 7^ Mj>^ because of their different evolution patterns. 

In addition, the D-tcrm contribution for the former causes a difference between the two 
diagonal entries. While the evolution of all other parameters in this scenario are practically 




(3) 



Csoft M?Jur\^ + {y.A^H^.Lu'i^ + h.c.) 



(4) 
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the same as in the MSSM, the right-chiral sneutrino mass parameter evolves [18] at the 
one-loop level as: 

Clearly, the extremely small Yukawa couplings cause M^^ to remain nearly frozen at the 
value mo, whereas the other sfcrmion masses are enhanced at the clcctroweak scale. Thus, 
for a wide range of values of the gaugino mass, one naturally has sneutrino LSP, which, for 
every family, is dominated by the right-chiral state: 

z/i = sin 9 + ur cos 9 (7) 

The mixing angle 9 is given as 

tan2g^ '^-"^^"f ^"^^,^-^-1 (8) 

which is suppressed by y^^ especially if the neutrinos have Dirac masses only. It is to be 
noted that all three (dominantly) right sneutrinos have a similar fate here, and one has near- 
degeneracy of three such LSPs. However, of the three charged slepton families, the amount 
of left-right mixing is always the largest in the third (being, of course, more pronounced for 
large tan/?), and the lighter stau (fi) often turns out to be the NLSP in such a scenario. ^ 

Thus the mSUGRA parameter set (mo, M1/2, Aq, sign{fi) and tan (3) in an 7?.-parity con- 
serving scenario can eminently lead to a spectrum where all three generations of right- 
sneutrinos will be either stable or metastable but very long-lived, and can lead to different 
decay chains of supersymmetric particles, as compared to those with a neutralino LSP. How- 
ever, as we shall see below, the deciding factor is the lighter sneutrino mass eigenstate of the 
third family, so long as the state fi is the lightest among the charged sleptons. 

All superparticles will have to decay into the lighter sneutrino of a particular family via 
either gauge interactions (such as tl — > Wu'^i) or Yukawa coupling (such as II — > H~Di or 
z/2 — > hPui). In the former case, the decay depends entirely on the Dl content of i>i, which 
again depends on the neutrino Yukawa coupling. The same parameter explicitly controls the 
decay in the latter case, too. Therefore, while the lighter sneutrinos of the first two famihes 
can in principle be produced from decays of the corresponding charged sleptons, such decays 



^We have neglected inter- family mixing in the sneutrino sector in this study. While near-degenerate 
physical states makes such mixing likely, the degree of such mixing is model-dependent, and does not generally 
affect the fact that all cascades culminate in the lighter stau, so long as the latter is the NLSP, which is the 
scenario studied here. 
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will be always suppressed compared to even three-body decays such as ei(/ii) — > e(/x)ffi 
(when the sleptons are hghter than all neutrahnos). For the f-NLSP, however, the only 
available decay channel is fi — W{H~)i>i, with either real or virtual charged Higgs. Both 
of these decay channels are driven by the extremely small neutrino Yukawa coupling. 

This causes the NLSP to be a long-lived particle and opens up a whole set of new 
possibilities for collider signatures for such long-lived particles, while retaining contributions 
to dark matter from the sneutrino LSP. The NLSP appears stable in collider detectors and 
gives highly ionizing charged tracks. 

Apart from a stau, the NLSP could be a chargino, a stop or a sbottom. The former 
is in general difficult to achieve in a scenario where the chargino and neutralino masses are 
determined by the same set of electroweak gaugino and Higgsino masses. The phenomenology 
of the long-lived stop NLSP [19], the likelihood of the corresponding signals being available 
at the early phase of the LHC, and the potential for the reconstruction of gluino masses have 
been discussed in an earlier work [12]. The stau NLSP, as we shall see below, offers a new 
opportunity to reconstruct both the lightest and second hghtest neutralino masses. 

It may be noted here that the region of the mSUGRA parameter space where we work is 
consistent with all experimental bounds, including both collider and low-energy constraints 
(such as the LEP and Tevatron constraints on the masses of Higgs, gluinos, charginos and 
so on as well as those from b — > 57, correction to the p-parameter, (g^^ — 2) etc.). Our 
choice of parameters in the mo — M1/2 plane would correspond to a stau-LSP without the 
right-sneutrino in the (super) particle spectrum. Such a situation would have been ruled 
out, had not the existence of the right-chiral neutrino superfield, with the right sneutrino 
at the bottom of the spectrum, been assumed [20] . However with the right-sneutrino as the 
LSP, we find this choice to be a preferable and well-motivated option. The i>i-LSP arising 
out of such a choice becomes a viable dark matter candidate, though not necessarily the only 
one. Using the formulae given by Moroi et al. in reference [8], the contribution to the relic 
density {ilh'^) is found to be about one order of magnitude below the acceptable value [9]. 
While this leaves room for additional sources of dark matter, the scenario presented here is 
consistent from the viewpoint of over-closure of the Universe. 

We focus on both the regions where (a) rrif^ > my-^ +mw, and (b) the above inequality is 
not satisfied. In the first case, the dominant decay mode is the two-body decay of the NLSP, 
fi i^iW, and, in the second, the decay takes place via a virtual W. However, the decay 
takes place outside the detector in both cases. Decays into a charged Higgs constitutes a 
subdominant channel for the lighter stau. 
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Furthermore, we try to identify regions of the parameter space, where neutrahnos de- 
caying into a tau and a stau can be reconstructed, through the reconstruction of the tau 
and the detection of the stau in the muon chamber. The rates for electroweak production of 
neutrahnos are generally rather low for this process. Therefore, the procedure works better 
when neutrahnos are produced from the cascade decays of squarks and gluinos. This is in 
spite of the additional number of jets in such processes, which may fake the tau in certain 
cases and comphcate the analysis of the signals. We are thus limited to those regions of 
the parameter space, where the gluino and squark production rates are appreciable, and 
therefore the value of M1/2 is not too high. 

With all the above considerations in mind, we concentrate on the lighter stau (fi) to be the 
NLSP with lifetime large enough to penetrate collider detectors like the muons themselves. 
Using the spectrum generator of ISAJET 7.78 [21], we find that a large mSUGRA parameter 
space can realize this scenario of a right-sneutrino LSP and stau NLSP, provided that rriQ < 
mi/2 and one has tanjS in the range ^ 25, the latter condition being responsible for a larger 
left-right off-diagonal term in the stau mass matrix (and thus one smaller eigenvalue). In 
Table 1 we identify a few benchmark points, all within a SUGRA scenario with universal 
scalar and gaugino masses, characterised by long-lived staus at the LHC. 

In the next section we use these benchmark points to discuss the signatures of the stau 
NLSP at the LHC and look for the final states in which it is possible to reconstruct the 
neutrahnos. 

3 Signal and backgrounds 

The signal which we have studied as a signature of stau NLSP and motivated by the possible 
reconstruction of the neutrahnos from the final state, is given by 

• 2tj + 2f {charged — track) + f]T + X 

where tj denotes a jet out of a one-prong decay of the tau, stands for missing transverse 
energy and all accompanying hard jets arising from cascades are included in X. 

We have simulated pp collisions with a centre-of-mass energy Ecm — 14:TeV. The pre- 
diction of events assumes an integrated luminosity of 300 fb"^. The signal and various 
backgrounds are calculated using PYTHIA (version-6.4.16) [22]. Numerical values of various 
parameters, used in our calculation, are as follows [7]: 
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Table 1: Proposed benchmark points (BP) for the study of the stau-NLSP scenario in the SUGRA 
with right- sneutrino LSP. The value of mQ and M1/2 (^re given in GeV. We have also set Aq = 
100 GeV and sgn{ii) = + for benchmark points under study. 



Mz = 91.2 GeV Mw = 80.4 GeV Mt = 171.4 GeV Mh = 120 GeV 
a-^{Mz) = 127.9 a,(Mz) = 0.118 

We have worked with the CTEQ5L parton distribution function [23]. The factorisation 
and renormahsation scale are set at /xf = /Ur = n^average- order to make our estimate 
conservative, the signal rates have not been multiplied by any K-factor [24], while the main 
background, namely, that from tt production, has been multiplied by a K-factor of 1.8 
[25]. The effects of initial state radiation (ISR) and final state radiation (FSR) have been 
considered in our study. 
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3.1 Signal subprocesses: 

We have studied all SUSY subprocess leading to the desired final states. Neutralinos are 
mostly produced in cascade decays of strongly interacting sparticles. The dominant contri- 
butions thus come from 

• gluino pair production: pp — > gg 

• squark pair production: pp — > qiqj,qiqi 

• associated squeirk-gluino production: pp — > qg 



Subprocesses 


BP-1 


BP-2 


BP-3 
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BP-5 


BP-6 


All SUSY 
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4143 
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28864 


15439 




1616 


3897 


11061 


19765 


27785 


14426 



Table 2: The number of 2Tj + 2f (charged-track)+^T +X events, satisfying our basic cuts, for 
J Ldt = 300 GeV, from various channel. 

In addition, electroweak pair-production of neutralinos can also contribute to the signal 
we are looking for. The rates are, however, much smaller [see Table 2]. Moreover, the 
relatively small masses of the lightest and the second lightest neutralinos (as compared to 
the squarks and the gluino) cause the signal from such subprocesses to be drastically reduced 
by the cut employed by us on the scalar sum of all visible prs. For example, for benchmark 
point 1 (BPl), one has less than 5% of the total contribution from electroweak processes. 

The production of squarks and gluinos have potentially large cross sections at the LHC. 
For all our benchmark points hsted in Table 1, the gluino is heavier than the squarks. Thus 
its dominant decay is into a squark and a quark. Xi being mostly B dominated the main 
contribution to Xi production comes from the decay of right handed squarks {(Jr) and its 
decay branching ratio into the f-r pair is almost 100% when it is just above the lighter stau 
in the spectrum. On the other hand, the X2 mostly 1^3-dominated, and therefore the main 
source of its production is cascade decay of left-chiral squarks (ql) ■ Such a X2 can also decay 
into a f-T pair. 

If one can obtain complete information on the four-momentum of the f and the r, it 
is thus possible to reconstruct both Xi and X2 using the final state mentioned above. The 

other two heavier neutralinos {x^ and xT) ? due to their low production rate and small decay 
branching ratios into the f-r pair, are relatively difficult to reconstruct. 
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3.2 Background subprocesses: 

The standard model background to 2Tj + 2f + comes mainly from the following sub- 
processes : 

• tt production, tt — > bWbW: 

Where two of the resulting jets can be faked as a tau-jet, and muons can come from 
the W's. One can also have a situation in which any (or both) of the b-quark decays 
semileptonically {b — > c/jlUh) and any (or both) of the W decays to a r — i/,- pair. Though 
the efficiency of a non-tau jet being identified as a narrow tau-like jet is small (as will 
be discussed in a later section), and it is very unlikely to have isolated muons from 
semileptonic decays of the b, the overwhelmingly large number of tt events produced 
at the LHC makes this subprocess quite a serious source of backgrounds. 

• Z°-pair production, ZZ — > 2t + 2//: 

This subprocess also gives an additional contribution to the background, when one Z 
decays into a tt pair and the other one into a pair of muons. 

• Associated ZH-production, ZH ^ 2r + 2fi: 

This subprocess, though have a small cross-section, can contribute to the background 
through the decay of if rr and Z — > ii/i, which can fake our signal as well. 

The additional backgrounds from ttW, ttZ and Z-l-jets can be suppressed by the same 
cuts as those described below. Also a higher non-tau-jet rejection factor and Z invariant 
mass cut can reduce the ttZ and Z-l-jets backgrounds considerably. 

3.3 Event selection criteria : 

In selecting the candidate events selected for neutralino reconstruction, we choose the two 
highest-pT isolated charged tracks showing up in the muon chamber, both with pt > 100 GeV, 
as stable staus. (See the detailed discussion later in this section.) The isolation criteria for 
the tracks are shown in Table 3. In addition to the pr cut, the scalar sum of all jets and 
lepton in each event is required to be greater than 1 TeV. It is clear from Figures 1 and 2 
that the standard model backgrounds are effectively eliminated through the above criteria. 
In addition, we require two r-jets with px > 50 GeV, and > 40 GeV for each event. The 
justification for both of these cuts is provided when we discuss r-tagging and reconstruction. 

The identification and 3-momentum reconstruction of the charged track at the muon 
chamber is done following the same criteria and procedure as those for muons. 
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Figure 1: pt distribution (normalised to unity) of the muon like track for the signal and the 
background, for all benchmark points. The vertical lines indicate the effects of a pt cut at 
100 GeV. 



In order to obtain the invariant mass of a tau-stau pair, one needs to extract informa- 
tion on the mass of the stable charged particle (the stau in our context). While standard 
techniques such as time delay measurement or the degree of ionisation produced has been 
suggested in a number of earlier works [26], we extract mass information from an event-by- 
event analysis which is reported in the next section. The efficiency for the reconstruction of 
staus has been taken to be the same as that of muons with px > 10 GeV in the pseudora- 
pidity range |?7| < 2.5, and is set at 90% following [27]. 
T-jet tagging and r-reconstruction: 

r-jet identification and r-reconstruction are necessary for both background reduction and 
mass reconstruction of the neutralinos. We have concentrated on hadronic decays of the r in 
the one-prong channel. ^. These are jet-like clusters in the calorimeter containing a relatively 
small number of charged and neutral hadrons. A r decays hadronically about 65% of the 
time, producing a r-jet. The momentum of such a jet in the plane transverse to the parent 
r is small compared to the r-energy, so long as the pt of the r-jet is large compared to the 
tau mass. In this limit, hadronic r decays produce narrow jets collinear with the parent r. 
The neutrinos that carry missing Et also have the same direction in this limit. This gives 

^We have not considered the leptonic decay of tau, as it is difficult to identify lepton coming from tau 
decay to that coming from cascade decay of other objects hke heavy quarks or W's 
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Scalar-sum of pt in GeV Scalar-sum of pt In GeV 

Figure 2: S|j9t| distribution (normalised to unity) for the signal and the background, for all 
benchmark points. 



one a handle in reconstructing the r's, if one selects events where no other invisible particle 
is likely to be produced. 

Given the masses of the SUSY particles in our benchmark scenarios, the r's produced 
out of neutralino decay are hard enough, so that one can simulate r-decays in the coUinear 
approximation described above. A detailed discussion on the procedure followed for com- 
plete reconstruction of a pair of r's is found in [28]. We have selected hadronic jets with 
Ex > 50 GeV as candidate products of r-decay. A rather conservative non-tau jet rejection 
factor of 20 has been assumed, while the identification efficiency of a true tau-jet has been 
assumed to be 50% following [29, 30]. 

To describe the procedure in brief, one can fully reconstruct the r by knowing Xr,^- (i = 
1,2), the fractions of the parent r-energy carried by each product jet. The two unknowns can 
be solved from the two components of the missing transverse momentum (Jjt) of a particular 
event. 

If pi^-,p'^i are, respectively the components of four-momentum of the parent r and the 
coUinear jet produced from it (i = 1,2), then 

P'H^ = Xr,^ P'r. (/i = 0, 1, . . , 4) (9) 

(as E't- ~ Ip^I, in the limit — ^ 0) 
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and one can write 



}t^{^- 1) Phi + - 1) Ph2 (10) 

This yields two conditions for x^-^^. Solving them, one obtains the r four- momenta as 

— * 

Phi/xrf^^. In practice, as will be discussed below, the recorded missing momentum, i^^?*^, is 
different from the true one, namely, This error can lead to unphysical solutions for 

the reconstructed r-momenta in some cases. Such a situation often arises when the two 
taus are produced back-to-back. This in turn means that the r-neutrino's are also produced 

— * 

back-to-back in the coUinear approximation. This reduces the magnitude of ^t, when errors 
due to mismeasurements can lead to unphysical solutions. This is sometimes avoided by 
leaving out back-to-back orientation of the two r-jet candidates, with some tolerance. In our 
analysis, a minimum value for {pd 40 GeV) and positivity of a^^r^/s have been imposed 
as necessary conditions, in order to minimise the number of unphysical solutions. Besides, 
Pt > 50 GeV for each r-jets ensures the validity of the coUinear approximation in r-decays. 
The T-identification efficiency and the jet rejection factor are also better optimized with this 
PT-cut [29, 30]. 

Of course, with a jet rejection factor of 1/20, one cannot rule out the possibility of QCD 
jets masquerading as r's, in view of the huge number of QCD events at the LHC. Such 
fakes constitute irreducible backgrounds to the r-reconstruction procedure. However, as we 
shall see in the numerical results presented in the next section, triggering on the rather 
strikingly spectacular properties of the quasi-stable stau-pair enables one to filter out the 
genuine events in the majority of cases. 
Reconstruction of ; 

It is evident from the above observations that the reconstruction of is very crucial for 
our study. The reconstructed differs considerably from true due to several reasons. 
The true ^^""^ is related to the experimentally reconstructed '^'^'^ by the following relation 

f^'' = + ^i;"'-"' + ^<o-5 (11) 

where corresponds to the total transverse momentum carried by the particles escaping 
detection in the range I77I > 5 and corresponds to the total transverse momentum 

carried by the particles in the range I77I < 5 with pt < 0.5 GeV^, which contributes to the 

— * 

true -^T- In addition to this, mismeasurement of the transverse momenta for jets, leptons 
3The threshold is 0.5 GeV for CMS and 1 GeV for ATLAS. 



14 



etc. alters the true by an sizable amount. This is due to the finite resolution of detectors, 
and is handled in theoretical predictions by smearing the energy /momentum of a particle 
through a Gaussian function. 

In our study, we have tried to reconstruct taking into account all the above issues. 
The missing transverse momentum in any event is defined as 

where the Ep^***^'^ consists of isolated leptons/photons/jets and also those objects which do 
not belong to any of these components but are detected at the detector, constituting the 
so called 'soft part' or the 'unclustcrcd component' of the visible momentum. We describe 
below the various components of the visible pr, and their respective resolution functions. 
Resolution effects: 

Among the finite resolution effects of the detector, taken into account in our analysis, 
most important are the finite resolutions of the'elcctromagnetic and hadron calorimeters, 
and the muon track resolution. Since the kind of final state we are dealing with does not 
require any isolated electrons/photons, we have not considered electron or photon resolution. 
The electrons/photons which are not isolated but have Et > 10 GeV and |?7| < 5 have been 
considered as jets and their resolution has been parametrised according to that of jets. Jets 
have been defined within a cone of Ai? = 0.4 and > 20GeV using the PYCELL fixed 
cone jet formation algorithm in PYTHIA. Since the staus are long hved and live a charged 
track in the muon chamber, their smearing criteria have been taken to be the same as those 
of isolated muons. Though one can describe the resolution of the track of staus and muons 
by different resolution functions (as nif » m^), one does not envision any significant devi- 
ation in the prediction of events via such difference. Therefore, in the absence of any clear 
guidelines, we have treated them on equal footing, as far as the Gaussian smearing function 
is concerned. The tracks which shows up in the muon chamber, but are not isolated, having 
Et > lOGeV and \r]\ < 2.5, have been considered as jets and smeared accordingly. All the 
particles (electron, photon, muon, and stau) with 0.5 < Et < lOGeV and \ri\ < 5 (for muon 
or muon-like tracks, I77I < 2.5), or hadrons with 0.5 < Et < 20GeV and I77I < 5, which consti- 
tute 'hits' in the detector, are considered as soft or unclustered components. Their resolution 
function have been considered separately. Wc present below the different parametrisation of 
the different component of the final state, assuming the smearing to be Gaussian in nature. ^ 



Although a "double Gaussian" smearing is followed in more realistic detector simulations, our illustrative 
study is now substantially affected by such considerations. 
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• Jet energy resolution 



a{E) /E = a/\fE ®h®c/E (12) 

where 

a= 0.7 [GeVV2], b= 0.08 & c= 0.009 [GeV] for I77I < 1.5 
= 1 = 0.1 = 0.009 1.5 < I77I < 5 

• Muon/Stau pr resolution : 



(^{Pt)/pt = a if pr < C (13) 

= a + b\og{pT/0 ifpT>C (14) 



where 

a= 0.008, b= 0.037 & ^ = 100 [GeV] for I77I < 1.5 
= 0.02 = 0.05 ^ = 100 1.5 < I77I < 2.5 

Soft component resolution : 



aiEr) = a^/E^E^' (15) 

with a ^ 0.55. One should keep in mind that the x and y component of E^^^ need to 
be smeared independently and by the same quantity. 

It is of great importance to ensure that the stable f leaving a track in the muon chamber is 
not faked by an actual muon arising from a standard model process. As has been mentioned 
in section 1, we have found certain kinematic prescriptions to be effective as well as reliable 
in this respect. In order to see this clearly, we present the pydistributions of the harder 
muon and the f-track in Figure 1. The f-pr clearly shows a harder distribution, owing to the 
fact that the stau takes away the lion's share of the pr possessed by the parent neutralino. 
Another useful discriminator is the scalar sum of transverse momenta of all detected particles 
(jets, leptons and unclustered components). The distribution in defined in the above 

manner, displays a marked distinction for the signal events, as shown in Figure 2. The cuts 
chosen in Table 3 have been guided by both of the above considerations. They have been 
applied for all the benchmark points, as also for the background calculation. 
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Cuts 


Basic Cuts 


^lep,stau ^ 
^hardest-jet ^ -,5 
^other-jets ^ g^y 

40 GeV <pT < 150 GeV 
\ri\ < 2.5 for Leptons, Jets & Stau 
ARii > 0.2, ARij > 0.4 
ARfi > 0.2, ARfj > 0.4 
ARjj > 0.7 


Cuts for 

Background 

Elimination 


E|pt| > 1 TeV 


Invariant mass difference 
of two nearby pairs 


|M?:"''-M?:"'^| <50 GeV 

1 TT TT 1 



Table 3: Cuts applied for event selection, background elimination and neutralino reconstruction. 

4 Numerical results and neutralino reconstruction: 

4.1 The reconstruction strategy 

Having obtained the r four-momenta, the neutrahnos can be reconstructed, once we obtain 
the energy of the f's whose three-momenta arc aheady known from the curvature of the 
tracks in the muon chamber. For this, one needs to know the f-mass. In addition, the re- 
quirements arc, of course, sufficient statistics, minimisation of errors due to QCD jets faking 
the r's, and the suppression of combinatorial backgrounds. For the first of these, we have 
presented our numerical results uniformly for an integrated luminosity of 300 fb~^, although 
some of our benchmark points requires much less luminosity for effective reconstruction. We 
have already remarked on the possibility of reducing the faking of r's. As we shall show be- 
low, a systematic procedure can also be adopted for minimising combinatorial backgrounds 
to the reconstruction of neutralinos. The primary step in this is to combine each such r with 
one of the two hardest tracks in the muon chamber, which satisfy the cuts listed in Table 
3. A particular r is combined with a heavy track of opposite charge. However, since neu- 
tralinos are Majorana fermions, producing pairs of r+f^ and r^f"^ with equal probability, 
this is not enough to avoid the combinatorial backgrounds. Therefore, out of the two r's 
and two heavy tracks, we select those pairs which give the closer spaced invariant masses, 
with \M?f''^ - M?f'''^\ < 50 GeV. The number of signal and background events, after the 
successive apphcation of cuts are fisted in Table 4. 

This finally brings us to the all-important issue of knowing the stau-mass. The f-mass 
can be reconstructed from the information on time delay (At) between the inner tracker and 
the outer muon system and the measured three- momentum of the charged track [26]. The 
accuracy of this method depends on the accurate determination of At, which can be limited 
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CUTS 






SIGNAL 






BACKGROUND 




BPl 


BP2 


BP3 


BP4 


BP5 


BP6 


ti 


ZZ 


ZH 


total 


Basic Cuts 


1765 


4143 


11726 


20889 


28864 


15439 


4129 


45 


6 


4180 


+Pt Cut 


1588 


3631 


9471 


15526 


20282 


9920 


210 


3 


1 


214 


+S|Pt| Cut 


1442 


3076 


6777 


9538 


11266 


5724 


63 








63 




408 


887 


1622 


2004 


2244 


858 


6 








6 



Table 4: Number of signal and background events for 2rj + 2f ( charged-track) +f^T + X final state 
with an integrated luminosity of 300 fh~^, considering all SUSY processes. The standard model 
Higgs mass is taken to be 120 GeV in background calculation. 



when the particles are highly boosted. We have followed a somewhat different approach to 
find the actual mass of the particle associated with the charged track. We have found this 
method effective when both pairs of rf come from x?x5 or y^y^. 

Solving for the f-mass: The actual f-mass can be extracted by demanding that the 
invariant mass of the two correct fr-pairs are equal, which yields an equation involving one 
unknown, namely, mf : 

Y^m? + b^ip,^^^ _ ymfTIp^.E^^ = p^^.p;, - (16) 

where the variables have their usual meanings and are experimentally measurable event-by- 
event. (fi,2 here denote the hghter fs on two sides of the cascade, and not the two f mass 
eigenstates.) 

The right combination is assumed to be selected whenever the difference between |M|"*^^ — 
j^po-tr2 1 minimum and differ by not more than 50 GeV, as mentioned earlier. It should be 
noted that the unambiguous identification of the right rf-pairs which come from decays of 
two neutralinos {xiXi or X2X2) does not depend on the actual stau mass. Thus we can use 
a 'seed value' of the stau mass as input to the above equation, in identifying the right rf 
combinations. We have used a seed value of nif ~ 100 GeV (motivated by the LEP limit on 
nif). The SM background has already been suppressed by demanding the px of each charged 
track to be greater than 100 GeV, together with T,\pt\ > 1 TeV. 

Once the right pairs are chosen using the seed value of the f-mass, we need not use that 
value any more, and instead solve equation 16 which is quadratic in m?. We have kept 
only those events in which at least one positive solution for m? exists. When both roots of 
the equation are positive, the higher value is always found to be beyond the reach of the 
maximum centre-of-mass energy available for the process. Hence we have considered the 
solution corresponding to the lower value of the root. The distribution of the solutions thus 
obtained has a peak around the actual f-mass. The f-track four-momenta are completely 
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constructed, using this peak value as the actual mass of the f-NLSP (see Figure 3). The fact 
that these peaks faithfully yield the f-mass (see Table-1), makes it unnecessary to extract 
this mass from 'fits'. 

This sets the stage fully for the reconstruction of neutralinos, the results of which are 
shown in Figure 4. For BPl, BP2 and BPS one can see that there is only one peak which 
corresponds to the Xv This is because the X2 production rate in cascade is relatively small 
for these points. For BP4 and BPS, on the other hand, we have distinct peaks for both Xi 
and X2- At BP6, however, we only have the X2 peak. This is due to the small mass splitting 
between Xi (-^t? ^ ^^9 GeV) and f (Mf = 124 GeV), which softens the tau (jet) arising 
from its decay, preventing it from passing the requisite hardness cuts. On the whole, it is 
clear from Figure 4 (comparing the peaks with the input values of the neutralino masses) 
that, in spite of adulteration by QCD jets that fake the r's, our event selection criteria can 
lead to faithful reconstruction of neutralino masses. 

One may still like to know whether a neutralino reconstructed in this manner is the Xi 
or the when only one peak is visible. This requires further information on the SUSY 
spectrum. For example, the information on the gluino mass, extracted from the effective 
mass (defined as (S|pt| + ^t)) distribution, may enable one to distinguish between the Xi 
and the X2: o^ce gaugino mass unification at high scale is assumed. 

4.2 LHC reach in the mo — M1/2 -plane: 

We have also identified the region in the mo — M1/2 plane, where at least one of the two 
lightest neutrahnos can be reconstructed. For this, we have scanned over the region of the 
mo — Ml/2 plane using the spectrum generator SuSpect (v 2.34) [31] which leads to a f LSP 
in a usual mSUGRA scenario without the right handed sneutrino [32]. Results of this scan 
are shown in Figure 5. The coloured (shaded) areas are consistent with all the low energy 
constraints like b S7, Bs fJ'^fJ' , ^(S'm ~ 2) LEP limits on the low energy 

spectrum. The value of tanf3 has been fixed at 30, and Aq = 100 has been chosen. 

The regions where reconstruction is possible have been determined using the following 
criteria: 

• In the parameter space, we have not gone into regions where the gluino mass exceeds 
^ 2 TeV. 

• The number of events satisfying \MfT- — Mpg^kl < 0.1. M^eafe at an integrated luminosity 
of 300 fb'^ must be greater than a specific number in order that the peak is said to 
be reconstructed. One obtains the blue (dark shade) region if this number is set at 
100. If the peak can be constructed from more than 50 events, the additional region, 
marked in green (light shade), becomes allowed. 
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Figure 3: The r mass peak as obtained from eventwise reconstruction as described i 
text, for all the benchmark points at luminosity 300 fb"^. 
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iUf, in GeV il/f, in GeV 



Figure 4: Mfr distribution for all the benchmark points at luminosity 300 fh . BPl, BP2 
and BPS show only the x? peak. Both the x? cind X2 peaks are visible for BP4 and BPS, 
while BP6 displays only the X2 peak. 
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200 I ' ' ' ' < 1 

100 150 200 250 300 350 

mO 

Figure 5: The region in the mo — M1/2 plane, where it is possible to reconstruct at least 
one of the neutralinos at the LHC, with tanf3 = 30 and Aq = 100. In the blue (dark shade) 
region, at least 100 events are predicted in the vicinity of the peak. The additional available 
region where 50 events in the vicinity of the peak are assumed to suffice for reconstruction, is 
marked in green (light shade). The entire region above the dashed line indicates the scenario 
where one has a z//j-LSP and a f-NLSP. 



5 Summary and conclusions 

We have considered a SUGRA scenario, with universal scalar and gaugino masses, where a 
right-chiral neutrino superfield exists for each family. We have identified several benchmark 
points in the parameter space of such a theory, where a right-sneutrino is the LSP, and a 
f mass eigenstate is the NLSP. The f, stable on the distance scale of the detectors, leaves 
a charged track in the muon chamber, which is the characteristic feature of SUSY signals 
in this scenario. We use this feature to reconstruct neutralinos in the rf channel. For this, 
we use the coUinear approximation to obtain the four-momentum of the r, and suggest a 
number of event selection criteria to reduce backgrounds, including combinatorial ones. We 
also suggest that the f mass may be extracted by solving the equation encapsulating the 
equality of invariant masses of two rf pairs in each event. We find that at least one of the 
two lightest neutralinos can be thus reconstructed clearly over a rather large region in the 
rriQ — Ml/2 plane, following our specified criteria. 
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